# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
results <- as.tibble(m1$est.att) %>%
mutate(n = row_number()) %>%
mutate(ci_90_low = ATT - 1.645 * S.E.,
ci_90_high = ATT + 1.645 * S.E.,
n = n-10.5)
ggplot(results, aes(n, ATT)) +
geom_hline(yintercept = 0, lty = 2, color = "gray60") +
geom_vline(xintercept = 0, lty  =2 ) +
geom_linerange(data = results,
aes(x = n, y = ATT, ymin = ci_90_low, ymax = ci_90_high),
size = 0.75,
color = 'gray80',
position = position_dodge(width = 1/2)) +
geom_point(position = position_dodge(width = 1/2), size = 2) +
labs(y = "ATT",
x = "Time") +
theme_minimal() +
theme(axis.title.y=element_blank(), legend.position = "none")
ggsave('./graphs/china/results_cluster2.pdf', width = 7, height = 3)
plot(m1, type = 'ct') +
theme_minimal() +
theme(legend.position = 'bottom') +
labs(title = "",
x = 'Treatment Time')
ggsave('./graphs/china/ct_cluster2.pdf', width = 7, height = 3)
knitr::include_graphics(c('./graphs/china/pre_treatment2.pdf',
'./graphs/china/results_cluster2.pdf',
'./graphs/china/ct_cluster2.pdf'))
df <- data_china %>% drop_na(IdealPointDistance)
df <- df %>% drop_na(polity2_dist, gdp_cap_dist)
df <- df %>% group_by(ctry1) %>% mutate(n = n())
df_sub <- df %>% filter(cluster == 3)
table <- df_sub %>% select(ctry1) %>% distinct()
knitr::kable(table)
ggplot(df_sub %>% filter(year %in% c(2004:2013)), aes(year, IdealPointDistance, color = ctry1)) +
geom_line() +
theme_minimal() +
theme(legend.position = 'none') +
labs(x = 'Year')
ggsave('./graphs/china/pre_treatment3.pdf', width = 7, height = 3)
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
results <- as.tibble(m1$est.att) %>%
mutate(n = row_number()) %>%
mutate(ci_90_low = ATT - 1.645 * S.E.,
ci_90_high = ATT + 1.645 * S.E.,
n = n-10.5)
ggplot(results, aes(n, ATT)) +
geom_hline(yintercept = 0, lty = 2, color = "gray60") +
geom_vline(xintercept = 0, lty  =2 ) +
geom_linerange(data = results,
aes(x = n, y = ATT, ymin = ci_90_low, ymax = ci_90_high),
size = 0.75,
color = 'gray80',
position = position_dodge(width = 1/2)) +
geom_point(position = position_dodge(width = 1/2), size = 2) +
labs(y = "ATT",
x = "Time") +
theme_minimal() +
theme(axis.title.y=element_blank(), legend.position = "none")
ggsave('./graphs/china/results_cluster3.pdf', width = 7, height = 3)
plot(m1, type = 'ct') +
theme_minimal() +
theme(legend.position = 'bottom') +
labs(title = "",
x = 'Treatment Time')
ggsave('./graphs/china/ct_cluster3.pdf', width = 7, height = 3)
knitr::include_graphics(c('./graphs/china/pre_treatment3.pdf',
'./graphs/china/results_cluster3.pdf',
'./graphs/china/ct_cluster3.pdf'))
df <- data_china %>% drop_na(IdealPointDistance)
df <- df %>% drop_na(polity2_dist, gdp_cap_dist)
df <- df %>% group_by(ctry1) %>% mutate(n = n())
df_sub <- df %>% filter(cluster == 4)
table <- df_sub %>% select(ctry1) %>% distinct()
knitr::kable(table)
ggplot(df_sub %>% filter(year %in% c(2004:2013)), aes(year, IdealPointDistance, color = ctry1)) +
geom_line() +
theme_minimal() +
theme(legend.position = 'none') +
labs(x = 'Year')
ggsave('./graphs/china/pre_treatment4.pdf', width = 7, height = 3)
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
data_china <- read_csv2("../data/final_data/china_data_with_cluster.csv") %>% filter(year > 2003)
data_china
data_china %>% distinct(cluster) %>% arrange(cluster)
test <- read_csv2('./final_data/china_data_with_cluster.csv')
test %>% distinct(cluster) %>% arrange(cluster)
cluster_data <- read_csv2('./final_data/cluster.csv')
knitr::kable(
list(cluster_data[1:30,], cluster_data[31:60,]),
caption = 'Comminuty Membership I',
col.names = c("Country", 'Cluster', 'Ties'),
booktabs = TRUE) %>%
kable_styling(font_size = 8)
knitr::kable(
list(cluster_data[61:89,], cluster_data[89:117,]),
caption = 'Comminuty Membership II',
col.names = c("Country", 'Cluster', 'Ties'),
booktabs = TRUE) %>%
kable_styling(font_size = 8)
knitr::kable(
list(cluster_data[1:31,], cluster_data[32:61,]),
caption = 'Comminuty Membership I',
col.names = c("Country", 'Cluster', 'Ties'),
booktabs = TRUE) %>%
kable_styling(font_size = 8)
knitr::kable(
list(cluster_data[62:93,], cluster_data[94:125,]),
caption = 'Comminuty Membership II',
col.names = c("Country", 'Cluster', 'Ties'),
booktabs = TRUE) %>%
kable_styling(font_size = 8)
df_sub
data_china <- data_china %>% filter(ctry1 %in% europe$countryname)
data_china <- read_csv2("../data/final_data/china_data_with_cluster.csv") %>% filter(year > 2003)
countries <- gisco_get_countries(
year = "2016",
epsg = "4326",
resolution = "10"
)
countries
countries <- countries %>% left_join(continents, by = "ISO3_CODE")
continents <- gisco_countrycode
continents
countries <- countries %>% left_join(continents, by = "ISO3_CODE")
countries <- subset(countries, CNTR_NAME != "Antarctica")
countries_small <- select(countries, CNTR_ID, continent) %>%
rename(Iso = CNTR_ID)
countries_small$geometry <- NULL
countries_small$countryname <- countrycode::countrycode(countries_small$Iso, "iso2c", "country.name")
europe <- countries_small %>% filter(continent == 'Europe')
europe
countries_small <- select(countries, CNTR_ID, continent, eu) %>%
rename(Iso = CNTR_ID)
countries_small$geometry <- NULL
countries_small$countryname <- countrycode::countrycode(countries_small$Iso, "iso2c", "country.name")
europe <- countries_small %>% filter(continent == 'Europe')
europe
europe %>% filter(eu = T)
europe %>% filter(eu = TRUE)
europe %>% filter(eu = 'TRUE')
europe %>% filter(eu == 'TRUE')
europe <- countries_small %>% filter(eu == 'TRUE')
data_china <- data_china %>% filter(ctry1 %in% europe$countryname)
df_sub <- data_china %>% drop_na(IdealPointDistance)
df_sub <- df_sub %>% drop_na(polity2_dist, gdp_cap_dist)
ggplot(df_sub %>% filter(year %in% c(2004:2013)), aes(year, IdealPointDistance, color = ctry1)) +
geom_line() +
theme_minimal() +
theme(legend.position = 'none') +
labs(x = 'Year')
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2 + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
results <- as.tibble(m1$est.att) %>%
mutate(n = row_number()) %>%
mutate(ci_90_low = ATT - 1.645 * S.E.,
ci_90_high = ATT + 1.645 * S.E.,
n = n-10.5)
ggplot(results, aes(n, ATT)) +
geom_hline(yintercept = 0, lty = 2, color = "gray60") +
geom_vline(xintercept = 0, lty  =2 ) +
geom_linerange(data = results,
aes(x = n, y = ATT, ymin = ci_90_low, ymax = ci_90_high),
size = 0.75,
color = 'gray80',
position = position_dodge(width = 1/2)) +
geom_point(position = position_dodge(width = 1/2), size = 2) +
labs(y = "ATT",
x = "Time") +
theme_minimal() +
theme(axis.title.y=element_blank(), legend.position = "none")
plot(m1, type = 'ct') +
theme_minimal() +
theme(legend.position = 'bottom') +
labs(title = "",
x = 'Treatment Time')
df_sub <- df_sub %>% drop_na(polity2_dist, gdp_cap_dist) %>% filter(ctry1 != 'France')
ggplot(df_sub %>% filter(year %in% c(2004:2013)), aes(year, IdealPointDistance, color = ctry1)) +
geom_line() +
theme_minimal() +
theme(legend.position = 'none') +
labs(x = 'Year')
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2 + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2 + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
results <- as.tibble(m1$est.att) %>%
mutate(n = row_number()) %>%
mutate(ci_90_low = ATT - 1.645 * S.E.,
ci_90_high = ATT + 1.645 * S.E.,
n = n-10.5)
ggplot(results, aes(n, ATT)) +
geom_hline(yintercept = 0, lty = 2, color = "gray60") +
geom_vline(xintercept = 0, lty  =2 ) +
geom_linerange(data = results,
aes(x = n, y = ATT, ymin = ci_90_low, ymax = ci_90_high),
size = 0.75,
color = 'gray80',
position = position_dodge(width = 1/2)) +
geom_point(position = position_dodge(width = 1/2), size = 2) +
labs(y = "ATT",
x = "Time") +
theme_minimal() +
theme(axis.title.y=element_blank(), legend.position = "none")
plot(m1, type = 'ct') +
theme_minimal() +
theme(legend.position = 'bottom') +
labs(title = "",
x = 'Treatment Time')
results <- as.tibble(m1$est.att) %>%
mutate(n = row_number()) %>%
mutate(ci_90_low = ATT - 1.645 * S.E.,
ci_90_high = ATT + 1.645 * S.E.,
n = n-10.5)
ggplot(results, aes(n, ATT)) +
geom_hline(yintercept = 0, lty = 2, color = "gray60") +
geom_vline(xintercept = 0, lty  =2 ) +
geom_linerange(data = results,
aes(x = n, y = ATT, ymin = ci_90_low, ymax = ci_90_high),
size = 0.75,
color = 'gray80',
position = position_dodge(width = 1/2)) +
geom_point(position = position_dodge(width = 1/2), size = 2) +
labs(y = "ATT",
x = "Time") +
theme_minimal() +
theme(axis.title.y=element_blank(), legend.position = "none")
plot(m1, type = 'ct') +
theme_minimal() +
theme(legend.position = 'bottom') +
labs(title = "",
x = 'Treatment Time')
df_sub <- df_sub %>% drop_na(polity2_dist, gdp_cap_dist)
ggplot(df_sub %>% filter(year %in% c(2004:2013)), aes(year, IdealPointDistance, color = ctry1)) +
geom_line() +
theme_minimal() +
theme(legend.position = 'none') +
labs(x = 'Year')
m1 <- gsynth(IdealPointDistance ~ treated + ctry1_fariss + log(ctry1_gdp_cap) + ctry1_polity2 + ctry1_polity2,
data = df_sub,
index = c("ctry1", "year"),
se = T,
force = 'time',
# CV = TRUE,
r = c(0, 5),
inference = "parametric", nboots = 1000,
parallel = FALSE)
knitr::opts_chunk$set(echo = TRUE)
# Load packages
library(tidyverse)
library(broom)
library(kableExtra)
library(igraph)
# Load data
network_base <- read_csv2('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/network_base.csv')
# The function takes the network_base data calculates the 5 percentile distance
# and compares the distance of each country to the 5th percentile. If it is  above it drops the country
id_function <- function(country, year) {
network_base %>%
filter(year %in% seq(t_year-10:t_year)) %>%
group_by(ctry1, ctry2) %>%
summarize(mean_dist = mean(IdealPointDistance)) %>%
drop_na() %>%
group_by(ctry1) %>%
# Here we vary the cut-off point
# Top 5th percentile
mutate(quantile_10 = quantile(mean_dist, 0.05)) %>%
filter(mean_dist < quantile_10)
}
c <- network_base %>% filter(treated_country == 1) %>%
select(ctry1, year) %>%
filter(year == 2013) %>%
distinct()
countries <- c$ctry1
t_year <- unique(c$year)
# Container
datalist = list()
# For loop for data with ID
for (i in countries) {
# ... make some data
dat <- id_function(i)
#dat$i <- i  # maybe you want to keep track of which iteration produced it?
datalist[[i]] <- dat # add it to your list
}
data <- bind_rows(datalist) %>% distinct()
data_saved <- data
# Here we vary the 0.25 threshold
data <- data %>% filter(mean_dist < 0.25)
ggplot(data, aes(mean_dist)) +
geom_histogram(bins = 50, color = 'black', fill = 'gray80') +
theme_minimal() +
labs(x = 'Average Distance',
y = 'Count')
ggsave('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/analysis/analysis/graphs/network/distribution-of-distances.pdf')
test <- knitr::kable(data %>% filter(ctry1 == 'Hungary'),
caption = "Table 1: Most proximate countries to Hungary",
col.names = c("Country 1", " Country 2", "Average Distance", "Fifth percentile"),
booktabs = TRUE,
format = "latex") %>%
kable_styling(font_size = 8)
# Relocate variables first to have connected entities
data <- data %>% select(-quantile_10) %>% relocate(ctry1, ctry2, everything())
# Create the network object
network_full <- graph_from_data_frame(d = data, directed = F)
# Calculate the clusters
#cluster_full <- cluster_optimal(network_full, weights = NULL)
cluster_full <- network_full  %>% cluster_walktrap()
# Storing the cluster information in a dataframe
# We will use this later to specify the clusters across the years
cluster_data <- as_tibble(cbind(cluster_full$names, cluster_full$membership)) %>%
arrange(V2) %>%
rename('ctry1' = 'V1',
'cluster' = 'V2') %>%
mutate(cluster = as.numeric(cluster))
degree <- degree(network_full)
number_of_ties <- enframe(degree) %>%
unnest
# Here we vary the the minimum number of ties
cluster_data <- left_join(cluster_data, number_of_ties,
by = c('ctry1' = 'name')) %>%
filter(value > 3)
ggsave('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/analysis/analysis/graphs/network/distribution-of-distances.pdf')
#write_csv2(cluster_data, 'C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/cluster.csv')
cluster_data <- read_csv2('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/cluster.csv')
# Create the voting network figure
V(network_full)$frame.color <- "white"
E(network_full)$arrow.mode <- 0
#### Reproduces Figure 4 in the Manuscript (note that the plot may look slightly different as the same network can be arranged in many different ways, however, the voting communities remain the same irrespective of the configuration of the network) ####
# Create the network figure
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
#dev.off()
#knitr::include_graphics('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
knitr::opts_chunk$set(echo = TRUE)
# Load packages
library(tidyverse)
library(broom)
library(kableExtra)
library(igraph)
# Load data
network_base <- read_csv2('C:/Users/CSteinert\OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/network_base.csv')
# Load packages
library(tidyverse)
library(broom)
library(kableExtra)
library(igraph)
# Load data
network_base <- read_csv2('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/network_base.csv')
rm(list = ls())
knitr::opts_chunk$set(echo = TRUE)
# Load packages
library(tidyverse)
library(broom)
library(kableExtra)
library(igraph)
# Load data
network_base <- read_csv2('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/network_base.csv')
# The function takes the network_base data calculates the 5 percentile distance
# and compares the distance of each country to the 5th percentile. If it is  above it drops the country
id_function <- function(country, year) {
network_base %>%
filter(year %in% seq(t_year-10:t_year)) %>%
group_by(ctry1, ctry2) %>%
summarize(mean_dist = mean(IdealPointDistance)) %>%
drop_na() %>%
group_by(ctry1) %>%
# Here we vary the cut-off point
# Top 5th percentile
mutate(quantile_10 = quantile(mean_dist, 0.05)) %>%
filter(mean_dist < quantile_10)
}
c <- network_base %>% filter(treated_country == 1) %>%
select(ctry1, year) %>%
filter(year == 2013) %>%
distinct()
countries <- c$ctry1
t_year <- unique(c$year)
# Container
datalist = list()
# For loop for data with ID
for (i in countries) {
# ... make some data
dat <- id_function(i)
#dat$i <- i  # maybe you want to keep track of which iteration produced it?
datalist[[i]] <- dat # add it to your list
}
data <- bind_rows(datalist) %>% distinct()
data_saved <- data
# Here we vary the 0.25 threshold
data <- data %>% filter(mean_dist < 0.25)
ggplot(data, aes(mean_dist)) +
geom_histogram(bins = 50, color = 'black', fill = 'gray80') +
theme_minimal() +
labs(x = 'Average Distance',
y = 'Count')
ggsave('C:/Users/CSteinert\OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/analysis/analysis/graphs/network/distribution-of-distances.pdf')
test <- knitr::kable(data %>% filter(ctry1 == 'Hungary'),
caption = "Table 1: Most proximate countries to Hungary",
col.names = c("Country 1", " Country 2", "Average Distance", "Fifth percentile"),
booktabs = TRUE,
format = "latex") %>%
kable_styling(font_size = 8)
# Relocate variables first to have connected entities
data <- data %>% select(-quantile_10) %>% relocate(ctry1, ctry2, everything())
# Create the network object
network_full <- graph_from_data_frame(d = data, directed = F)
# Calculate the clusters
#cluster_full <- cluster_optimal(network_full, weights = NULL)
cluster_full <- network_full  %>% cluster_walktrap()
# Storing the cluster information in a dataframe
# We will use this later to specify the clusters across the years
cluster_data <- as_tibble(cbind(cluster_full$names, cluster_full$membership)) %>%
arrange(V2) %>%
rename('ctry1' = 'V1',
'cluster' = 'V2') %>%
mutate(cluster = as.numeric(cluster))
degree <- degree(network_full)
number_of_ties <- enframe(degree) %>%
unnest
# Here we vary the the minimum number of ties
cluster_data <- left_join(cluster_data, number_of_ties,
by = c('ctry1' = 'name')) %>%
filter(value > 3)
ggsave('./replication_files/analysis/analysis/graphs/network/distribution-of-distances.pdf')
cluster_data <- read_csv2('C:/Users/CSteinert/OneDrive - Universität St.Gallen/Dokumente/unga_bri_project/replication_files/data/data/final_data/cluster.csv')
V(network_full)$frame.color <- "white"
E(network_full)$arrow.mode <- 0
#### Reproduces Figure 4 in the Manuscript (note that the plot may look slightly different as the same network can be arranged in many different ways, however, the voting communities remain the same irrespective of the configuration of the network) ####
# Create the network figure
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
#dev.off()
# Create the network figure
set.seed(1234)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
# Create the network figure
set.seed(12345)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
# Create the network figure
set.seed(123456)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
# Create the network figure
set.seed(1234567)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
# Create the network figure
set.seed(123456789)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
# Create the network figure
set.seed(123456789)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.85)
# Create the network figure
set.seed(123456789)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 10, vertex.label.cex=0.85)
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.45)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.65)
# Create the network figure
set.seed(123456789)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.65)
# Create the network figure
set.seed(123456789)
#pdf('C:\Users\CSteinert\OneDrive - Universität St.Gallen\Dokumente\unga_bri_project\replication_files\analysis\analysis\graphs\network\full-network.pdf')
plot.igraph(network_full,  vertex.color=membership(cluster_full), vertex.size = 7, vertex.label.cex=0.85)
